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Sampling 


. V(x, y,z):x+y+z=0 (X1, Vi, Z1) and (x2, y2, Z2) are the 
= Coordinate system co-ordinates of the two hexagons. 


o hexagonal grid Drud ((X1: V1. X1), (X2, Y2, x2)] 


o SHCS: symmetrical hexagonal = Valle — x2)? + (y1 — yz + (a — 22)" 


coordinate system in (C) 


Darid{(x1. y1, x1). (x2, y2, x2)] = max(Ixi — xol. ly1 — y21, 121 — 221) 
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V(x,y,z):x+y+z=0 Qu, Y1, 21) and (x2, y2, 27) are the 
co-ordinates of the two hexagons. 


Drud ((X1. y1. x1). (32. y2. X2)] 
ES V Alte — x)? + (y1 —yay + (z1 — z2)”] 


Deridl(X1. y1, x1). (x2. y2, X2)] = max(|x1 — x2|, |y1 — yal, |Z1 — zal) 
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Sampling 


= Hierarchical grids 


o neighbours at finer scale become 
focal cell or centroids for coarser scale 


o smooth out or simplify some grids 


* dynamic grid resolution 


o 0,r can be used to find out current resolution scale 


Alternating CW, CCW 
19.1? rotations of 7 
children 1/7th the area 
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Sampling 


= Checkerboard effect 


o due to uniform non-optimal 
square grid sampling 
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0 otherwise 
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= Aliasing 
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" Aliasing Lin, y) = f(x, y) comb(z, y, Az, Ay) 


f(z,y) 


comb(z, y; Az, Ay) x Ka? 2:5 ó(z S mAz, y — nAy) 
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š(x,y) = 0 otherwise 
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sampling 


= Aliasing fs(z,y) = f(z, y) comb(z, y, Az, Ay) 


F (Wr, Wy) = Wr,Wq, Kä X F(w; — pus, Wy — QWu,) 
p q 


ad 0 wy 
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Fs (we, wy) F,(wz, wy) 
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Sampling 


Saraz(t, z) continuous bandlimited function 


JEŻ) 


= Sampling theorem 


o f(t,z) can be recovered fully 
with zero error from its samples 


o iff grid is ‘sufficiently’ dense 


o just change of variables to 
simplify notations 


oo 


SATAZ (t, z) = > O Kit — mAT, z == nAZ) 


m=-—00 n=—00 


F(u,v)=0 for |u| = puma and lvl = vma ... band limits 
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Sampling 


continuous bandlimited function 
= Sampling theorem 


o f (t, z) can be recovered fully 
with zero error from its samples 


Image «—————— f(t, z) 


o No info is lost in the image if it is obtained via sampling at rates 
greater than twice the max freq content of f (t, z) in both u, v directions. 


00 


Saraz (t, Z) = bu > d(t — mAT,z — nAZ) 


m=-—00 n=—00 


F(u,v) = 0 for p = Mmax and || = Vmax ... band limits 
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= Erroneous effects 
o square grid sampling 
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Sampling 


Erroneous effects 
square grid sampling 


input 8x8 sq grid 8x8 mean 
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Ouantization 


= Ouantizer 
o SISO - scalar quantizer 
o mappings [xi-,,xi) > yi 
o what are the unknowns? 


x € [tę, tk+1) => Q(z) — r 
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Quantizer 
1 c.ca 


Input signal: X 


Quantizer input 


Output levels 
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Ouantization 


= Uniform guantizers 


o all ranges divided egually 
with A = Itx, tx+1) intervals 


o deadzone 


Midtread guantizer Midrise guantizer 
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Ouantization 


= Non-uniform guantizers 


o ranges divided via predefined 
function which gives A intervals 


Logarithmic guantizer 
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Logarithmic quantizer for image filter weights 


linear quantizer log, quantizer 
3000 ~ TTT 3000 COE 
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= Non-uniform guantizers Original and guantized signal 
o can we design optimal guantizer? 
o optimal in the sense to minimize error (which error?) 
o input: x; = t; :thresholds 
o output: y; = 7; :reconstructions 


o signal distribution is known: p(x) 
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= MSE 


D= f (c-Q(2))?p(a)da 
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ky Tk k—1 “tk 
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E Tp" ie L= looses 


OD o L rou 
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Ouantization 


L tk+1 
minimize D= >| (x — rp)*p(x)dx 
subject to tı = a, trai =b, tk < fait, 
ip Ty" ie KS leone d£ 


Fundamental theorem of calculus: 
derivative with accumulation function 


aD ð L (tu 
Ot, > U (z — rj)^p(z)dx (c initial const., f (t) is cts in open interval) 
 3=1'7 


dt), 
A(t) = J fax 
A'(t) = f(t) 
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Ouantization 


L 


tk+1 
minimize D= >| (x — rp)*p(x)dx 


subject to tı =a, tr41= b, 
tk < Tk S tk+1 


SP d < tj+1 2 
Ot At); 2 JA (x — rj) p(z)dz 


= (ty — rk)” p(ti) + (ty — Th-1)? (th), k=2,..., 


— d ( [M 
Leibniz integral rule: dz ( 


a(x 
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tę Ui. 
ee ee 


Fundamental theorem of calculus: 
derivative with accumulation function 
(c initial const., f (t) is cts in open interval) 


A(t) = | fax 
A'(t) = f(t) 


z) d d b(z) ð 
| pea = fee) 506) ` (at) ale) + | | Ia 
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Ouantization 


L 


tk+1 
minimize D= >| (x — rp)*p(x)dx 


subject to tı = a, trai =b, tk < fait, 
ip Ty" ie KS leone d£ 


OD qe OD 
— = — d Së jtl 2 
Es TÈ 2, EDO = SH x VCI 


Ot}, 
= 2/7 "Ge —ry)p(z)dz, k = 1,...,L 
k 


E: - SCH + (tk — ry i 4p). k=2,...,1 


Fundamental theorem of calculus: 
derivative with accumulation function 
(c initial const., f (t) is cts in open interval) 


A(t) = JA f(x)dx 
A'(t) = f(t) 
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Ouantization 


L tk+1 
minimize D= >| (x — rę) p(x)dz 
subject to: tı =a, Tra =b, tk < fait, 
ip Ty" lie L= Leong hs 


OD Y [^ OD 
OE i © — r3)?p(z)dz sce i T — r;)?p(z)d 
Ory, E At), SH 3 (x Tj) p(z)dz 
= -2 igi "(x — rę)p(z)da, k=1,...,L = E: ECH + (tk — ry i 4p). k=2,...,1 
k 
Assuming p(x) > 0 for each z € [a, b]: Fundamental theorem of calculus: 
í derivative with accumulation function 
jg" zp(z)dz (c initial const., f(t) is cts in open interval) 
Tk “WEI Der è o = L,..., L 
Si, p(a)de A(t) = JA f(x)dx 
be = TR PT, NONE, A' (t) = f(t) 
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Ouantization 


= Lloyd-Max Quantizer poh zp(z)dz 
ak GEE be 
o optimal MSE quantizer E AE T C (1) 
i [ER pla)dz 
Tk-1 Tk _ HAT > .(2) 


2 2 i 
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Ouantization 


= Lloyd-Max Quantizer poh zp(z)dz 
as SESJĘ 
o optimal MSE quantizer KOSE TE pne „(1) 
Se, p(z)dz _ 
rr tro _ Tk-1tTk p O 


2 2 | 


a Pseudo-code 


: pick initial values for t (uniform grid) 
: find r values using (1) 

: find new t values using (2) 

: repeat till both t, r converge 
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Ouantization 


= Properties of optimal guantizer 


o E|Q(z)) = 35, mp _ fep ap(z)dz 


| LP) p(z)dz 
Elx] 


tk+1 
Pk = / p(x)dx 
` k 


o Bla —Q(x)|=0 
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Ouantization 


= Properties of optimal guantizer 


o  Elr-Q(rz)Q(x) = - Ef" (x — ry)rkp(x)da E Kaes 
E p(x)dx 
Pk = | m p(z)dz 
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Ouantization 


= Properties of optimal guantizer 


tk+1 tk+1 l I 
o  Elr-Q(x)Q(x) = X» (x — ry)rkp(z)da E ua 
Bon It plz)dz 
tk 
=0 = | TD o(z)dz 
Jt, 
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Ouantization 


= Properties of optimal guantizer 


tk+1 tk+1 l I 
o  Elr-Q(r)Q(x) = X» (x — ry)rkp(z)da E ua 
Bon It plz)dz 
tk 
=0 = | TD o(z)dz 
Jt, 


o error is uncorrelated with the quantizer's output 
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Quantization 


a Lloyd-Max example 


8 bpp 
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Quantization 


Lloyd-Max example 


8 bpp 6 bpp 
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Quantization 


Lloyd-Max example 


8 bpp 6 bpp 4 bpp 
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Conclusion 


- Sampling 


- Quantization 
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Sampling 
" Squares 
u Hexagonal 


" Aliasing 


Conclusion 


Quantization 


- Sampling = Uniform 
= Non-uniform 
- Quantization * Optimal 
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